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I. INTRODUCTION 



As it is well known, the spontaneous breakdown of the chiral symmetry in massless QCD 
is signaled by the appearance of a non-vanishing order parameter, the chiral condensate 
Xl =< >• This quantity is intrinsically non perturbative and its determination is not 
at all straightforward. Since numerical computations on the lattice yield -at a given fixed 
lattice spacing- only a bare result, one needs renormalization in order to use the result for 
phenomenological applications. The results of several recent studies 0, 0, O] evidence 
by now the dependence of the chiral condensate on both the scale and the renormalization 
scheme but allow for the computation of a scale independent Renormalization Group In- 
variant (RGI) chiral condensate jsj. Recently, an interesting analytical evaluation of the 
quark condensate in one-flavor massless QCD from the value of the gluino condensate in 
SUSY Yang-Mills theory by means of orientifold large-iV expansion has appeared in the 
literature 6|. 

In this paper we compute the bare chiral condensate in the Hamiltonian formulation of 
lattice QCD using staggered fermions in the strong coupling limit of the theory for any value 
of Nc; the lattice chiral condensate is evaluated using as the expansion parameter and, 
then, taking the 't Hooft limit (large Nc, with g'^Nc fixed). The extrapolation to the con- 
tinuum limit is carried out by means of Pade approximants. The tools used for our analysis 
were developed in jl| ; here we shall present only the result for the chiral condensate referring 
the reader to Q] for the details of the method. 

The result obtained for the chiral condensate xl is compared to the recent numerical deter- 
mination of the RGI chiral condensate pj]. In particular, we estimate the scale independent 
and dimensionless ratio between xl and the cube of the mass of the p meson, which has been 
obtained in a strong coupling large Nc calculation in 1], and compare it to the numerical 
value of the ratio between the scale independent RGI chiral condensate obtained in ^ and 
the experimental value of the cube of the mass of the p meson; the check allows for an 
estimate of the lattice light velocity t on the lattice, t turns out to be equal to 1.025, in 



t in turn confirms the accuracy 
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excellent agreement with the expected result of 1. This resu 
of our previous evaluation of the mass spectrum of mesons 

The strong coupling limit, where the hadrons are automatically confined, is a natural 
starting point for the study of the long-distance, non-perturbative features of the chiral 
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symmetry breaking in massless QCD. Even if, at weak coupling, the Wilson formalism yields 
an accurate analysis of the U{1) anomaly and PC AC one cannot use it to compute the 
chiral condensate at strong coupling since it explicitly breaks the chiral symmetry on the 
lattice. 

The staggered fermion formalism yields, instead, a remnant of chiral symmetry on the 
lattice, which is the invariance of the theory under translation by a single link; thus, a 
spontaneous breakdown of chiral symmetry may be evidenced also at strong coupling via 
a non vanishing chiral condensate. Moreover, the staggered fermion formalism is known to 
yield good results in the strong coupUng evaluation of the hadron spectrum |1 Q and 
of chiral condensate of lower dimensional models [10|, llH Other types of lattice fermions 
such as domain-wall or overlap fermions are expected to suffer both doubling and explicit 
breaking of chiral symmetryjl2?] at strong coupling. 

Str ong ly coupled lattice gauge theories are intimately related to quantum spin sys- 
tems in particular lattice gauge theories with staggered fermions exhibit interesting 
similarities with condensed matter systems. For example, it is well known that the quantum 
spin-1/2 Heisenberg antiferromagnet is eq uivalent to the strong coupling limit of either a 
f/(l) or SU{2) lattice gauge theory For the gauge group t/YA/'c), one may estab- 

lish the equivalence with a spin-A/c/S Heisenberg antiferromagnet 171. Since the seminal 
work of 't Hooft |l8l. Il9| the large Nc limit (with g'^Nc fixed), has played an increasingly 
rmportarrt role in studying gauge theories rn the contrnuunr. on the lattrce Q Q and the 
duality between gauge and string theories |21|,|22|. Furthermore^in a recent paper m, we 
have generalized the strong coupling calculation by Banks et al. jSi] to an arbitrary number 
of colors Nc and we have used the 't Hooft limit to investigate some features of the meson 
spectrum of strongly coupled lattice QCD. Our results imply that, also at strong coupling, 
the 't Hooft limit offers a very accurate method to consistently determine both the QCD 
spectrum and the chiral condensate. 



II. STRONG COUPLING EVALUATION OF THE QCD CHIRAL CONDENSATE 
IN THE 'T HOOFT LIMIT 

In the Hamiltonian formulation of lattice QCD with staggered fermions j^l time is a 
continuous variable and space is discretized on a 3-dimensional cubic lattice with M sites. 
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labeled by r = {x,y,z); with x, y and z integers. The lattice Hamiltonian with one flavor 
of massless quark may be written as the sum of three contributions 



H = He + Hg + Hm, (1) 



where 



He = l-^j:Er,nr (2) 

[r,h] 

H, = ^Y. Tl{n)^A{r + n)UAB[r^ n]^B(r) + h.c. = + Hi (3) 

[r,n] 

Hm = — TT' \Tr{U\r^n]U\f + n^m]U^\f + m^h]U^ (4) 

[r,ra,r?ij 

are the electric fleld Hamiltonian, the interaction Hamiltonian between quarks and gauge 
flelds and the magnetic Hamiltonian, respectively, t is the lattice light velocity. The sums 
Y^ir^h] ci-re extended to the lattice links, whereas Y.[r,h,Th] is a sum over the plaquettes. 
n = X, y, 5 is the unit vector in the n direction and 

//(x) = r^{y) = i-ir, V{^) = (-1)^ (5) 



are the Dirac a matrices for staggered fermions |23l]. The gauge fleld U[f,n\ is associated 
with the link [f,n\ ant it is a group element in the fundamental representation of SU{Nc)- 
In the strong coupling expansion the electric fleld Hamiltonian H^, is the unperturbed 
Hamiltonian while the interaction Hamiltonian between quarks and gauge flelds Hg, Q, 
and the magnetic Hamiltonian H^, Q, are treated as perturbations. An important feature 
of the Hamiltonian (Q) is its invariance under translation by a single link which plays the 
role of a discrete chiral symmetry [8|. It takes even sites into odd sites (even (odd) sites 
are those with x + y + z even (odd)). The pertinent lattice transformations have the form 
\['(r) \[^(r + x){-iy, \[^(r) \&(r + y){-iy and \I'(r) ^^(r + z){-lY. In momentum 
space the last equation can be written as g — > e*^"=75r3g, which in the continuum limit, where 
kz is inflnitesimal, becomes q — > 75X3^. The other two transformations yield q 75T2g, 

He has two degenerate ground states corresponding to those of the spin Nc/2 antifer- 
romagnetic Ising model. One state may be obtained from the other by interchanging odd 
and even sites; choosing one of these two vacua leads to the spontaneous breakdown of the 
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chiral symmetry. In the thermodynamic hmit the two ground states are not mixed to any 

I I 

finite order of perturbation theoryllO] and thus, in the perturbative expansion, one only 
has to consider diagonal matrix elements and, consequently, perturbation theory for non- 
degenerate states. The ground state energy has been evaluated in the strong coupling regime 
in m up to the fourth order in the strong coupling expansion and for a generic value of Nc- 
The order parameter evidencing discrete chiral symmetry breaking is the mass operator 
M. = il){r)il){r) since it acquires a nonzero expectation value giving rise to the chiral con- 
densate. In the staggered fermion formalism the pertinent lattice operator is given by 

^ = E(-l)^^'^^^i(^1^A(r). (6) 

r 

One has to evaluate the expectation value of M. on the perturbed states |po > generated by 
applying Hq to the ground state |0 >. One has |j9o >= |0 > +|po^^ > +bo^'' >, where 



The lattice chiral condensate is then given by 

_ < 0|A<|0 > + < p^^^l-Mb^^^ > + < vf\M\vf > 



where 

<,>= 



n/rff/[f,n]l(,) 



is the inner product in the full Hilbert space of the model. dU is the Haar measure on the 
gauge group manifold and (, ) the fermion Fock space inner product; IIq is the projection 
operator projecting onto states orthogonal to |0 >. 

Xl may be more conveniently computed by constructing an eigenstate of Hf. and using 
it to evaluate the functions f{He) appearing in Eq.ijH}. In order to show the method used 
to evaluate xl 3] we will concentrate on the denominator of Eq.ijH)). Since the vacuum 
state |0 > is a singlet of the electric field algebra, one has ii^"[r, fi,] |0 >= which, in turn, 
implies that He\0 >= 0. Using the left action of the Lie algebra generated by the electric 
field E"'[r,h] on U[r,h], He\0 >= and putting the commutator [ife, ^]] |0 > in place 
of HeU[r,n]\0 >, one finds 

H,U[f,n]\0 >= ^C2iNc)U[f,n\\0 >, (10) 



where C2 = {Nq — l)/2Nc is the Casimir operator of SU{Nc)- U[f,n\\0 > is then an 
eigenstate of He with eigenvalue g^C2{Nc)/2a. Consequently, 



< 0\H„ 
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HJO >-- 



8a^ 



< 



oihIhjo 



> . 



(11) 



Taking into account eq. pi|) after integration over the link variable U 2J| and the fact that 
< 0|0 >= 1, one finds 

(12) 



(2) (2) 

To compute < Po \pq > one needs to construct suitable eigenstates of the unperturbed 
Hamiltonian containing two and four link variables U[f,n\ Using the eigenvalues of 
the /(ife)'s on these eigenstates one finds 



4a^ 

7cl 



7cl 



< o\h,hIh,hI\o > + 



2a2 



< 0\HgH,HlHl\0 > 



OTSJ O 

+ ^ 9^2 < 0\H,H,Y.'^^AmUr,n]^B{r + n) 

y^^c - 2) ^^^^ 

X ^^(f)f/^o[f,n]^B(r + n)|0 >) 

+ ^ < 0| UAB[r,h]UBc[r + n,m]Ulij[r + m,h]ulj^[f, 

" 2 [r,h,'rfi] 

X E UefIt', ipFclr' + I, kpl^lr' + k, /]f/^s[f', A;]|0 > 

[r'lk] 



m\ 



(13) 



and the integration over the link variables leads to 



<P^o^\p^o^> 



+ 



2 



-5Nc - ^ 1 iV + '-^N'^ + 



Nl 2 , Nc{Nc-lf 
2 2{Nc-2Y 



g'cr 



(14) 
(15) 



The expectation values of the mass operator M. are computed using a similar procedure and 
they are given by 

Nc 



< Q\M\Q >-- 
<pi^^\M\pi'^ >- 



2a3 



g^CW 



N, 



-6Nc + —N 



, 17 , M 
60Nc + QNl - —NlN - -^N 



(16) 
(17) 



C at2 



N'c ,r2 , NciNc - If 



;-6 + ^-fN) 



{Nc - 2)2 



Nc N 



(18) 
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The last contributions in eq. ()15|) and in eq. ljlSp come from the magnetic term in \pq'^ >. 
Using eqs. p2llj!Hj) . one finds a non- vanishing chiral condensate for any finite Nc- For brevity, 
we do not write the expression of xl for a generic value of Nq, since it may be easily inferred 
from ip fT^ . It is easy to see that, as it should be, all the dependence on the number of 
lattice links disappears from the evaluation of xl up to the order t^/g^ and this is a very 
non trivial check of our calculations. 

Consider now the 't Hooft limit where q'^Nq is rescaled to g'^ and Nq is then sent to infinity. 
One gets for the bare lattice chiral condensate 

Nc 



Xl = o 



1-24:6 + 864e^ 



(19) 



where e = t^ / g^. 



III. COMPARISON WITH THE RGI CHIRAL CONDENSATE 

To compare the results of the strong coupling expansion with the numerical value of the 
scale independent RGI chiral condensate determined in P] one needs to construct a scale 
independent quantity. This can be done by considering for example the ratio between xl 
which has dimension (mass)^ and the mass cube of a meson. In a previous work ^\ we 
derived, in strong coupling, the series expansions for the masses of the low-lying states in 
the meson spectrum obtaining a good agreement between the experimental values for the 
meson mass ratios and our lattice results. For the p meson we found 



nip 



9_ 

a 



- + 6e - 203e2 



(20) 



.4 

In order to eliminate the dependence on a^, one may then consider the ratio between x = 
-Xl/Nc and 

•A/ 1 

2 



m^p g^ (i + 6e - 203e^ 
To compare this strong coupling result with the numerical value of the RGI chiral condensate, 
the series derived in the strong coupling regime (e = t'^/g^ ^ 1), need to be extrapolated to 
the region in which e ^ 1, which corresponds to weak coupling. To make this extrapolation 
possible it is customary to use the Fade approximant method which allows one to extrapolate 
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a series beyond its convergence radius. Using = e^/^/t^, from (j2ip one gets 

4/3 



mi 



4/3 



e 



i - 24e + 864e2 
i + 6e - 203e2 



(22) 



Using the [0, 2] Pade approximant, one has 



A. 

mi 



4/3 



For 



ooj, eq. 



becomes 



X 



m'i 



4/3 



2^/3 X 64 
F 1 + 160e + 7632e2 ' 



' 2^/3 X 64 \ l_ 
7632 / t^- 



(23) 



If one takes for the chiral condensate x the RGI numerical value obtained in and 
value of the mass of the p meson its experimental value, equating their ratio to 
obtains 



^ 2^/3 X 64 \ l_ 
7632 j 



(Xnum ± Ax num ) 

((0.771 ±0.0009))= 



I 4/3 



where Xnum is the numerical value of the chiral condensate of ref. 

Xnum ± Ax = (0.0167 ± 0.0031)Gev^ 
^From eq. (j25j) one gets for t the following value 

t = 1.025 ±0.177 



(24) 

for the 

, one 

(25) 
(26) 

(27) 



The result obtained is close to the expected value of 1, within the error, thus showing 
very good agreement between our results for the ratio Xi/'^p ^^id the one obtained from 
the results of 0. In the r.h.s. of eq.(j2SI) the error on the numerical value of xl is due to 
statistical effects while the error on the mass of p is experimental. Of course, for lattice 
QCD we can only check the consistency of the strong coupling determination of the chiral 
condensate by comparison of our result with the numerical value of the RGI chiral condensate 
determined in A comparison with the "true value" of the continuum chiral condensate is - 
for QCD- not only impossible but also senseless due to the manifest scale and renormalization 
scheme dependence of the lattice determinations of this quantity . As we shall see 

immediately after, in the simpler case of the Schwinger model, the strong coupling evaluation 
of the scale independent ratio between the chiral condensate and the pseudoscalar excitation 
mass matches very well the exact continuum value. 
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IV. CONCLUDING REMARKS 



The reliability of the method used to estimate the accuracy of our computation of the 
chiral condensate in QCD may be usefully tested in the one flavor Schwinger model for 
which the value of the chiral condensate and the mass of the pseudoscalar boson in the 
continuum are known exactly. In fact, one may compute the ratio between the lattice chiral 
condensate and the lattice pseudoscalar mass and then equate it to its continuum exact 

n 

counterpart. For the sake of clarity, we briefly recall the results of Berruto et al. for the 
evaluation of the lattice chiral condensate and the pseudoscalar mass of the strongly coupled 
one-flavor Schwinger model in the Hamiltonian approach with staggered fermions. In the 
continuum gauge model the chiral symmetry is broken by the anomaly. In the Hamiltonian 
lattice formulation of the Schwinger model the axial symmetry is spontaneously broken via 
a non zero expectation value of the chiral condensate. In the continuum theory it is well 
known that [25] < ■?/'(a;)^/'(x) >= —e'^Cc/ (271-^7?), where 7 = 0.577... is the Euler constant 
and m = ed \^ is the mass of the pseudoscalar excitation. 

A lattice Hamiltonian which, in the continuum limit, reduces to the Schwinger Hamiltonian 
is 

Hs=^^T.El-^ T.i^l+ie'^^^. - V^ie-^^^.+i) ^ + H,, (28) 

where the fermion fields are defined on the sites, x = — y,— y + gauge and electric 

fields. Ax and E^, on the links [x,x + 1]; is an even integer. The coefficient t of the 
hopping term in (j^H|) plays the role of the lattice light speed. In the naive continuum limit, 
ei = Cc and t = 1. 

In the strong coupling limit the electric field Hamiltonian is the unperturbed Hamiltonian 
while the hopping Hamiltonian Hp is treated as a perturbation. As in 3+1 dimensional QCD, 
there are two degenerate gauge invariant ground states which, in the Coulomb gauge, have 
the form 

i^>= n ^iio>> \x>= n ^iio>- (29) 

x=even x=odd 

To the fourth order in a = t/2e\a?, the perturbative expansion for the ground state energy 
is given by [l^ 



= Ef + a^Ef + a^Ef = ^ - ANa" + 192iVa^ (30) 
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On the lattice the pseudoscalar boson of the continuum Schwinger model is provided by 
the operator \6 >= I / {VN) J2x=ii'^x^^'^'^x+'i- ~^ i'x+i^'^^i'x)- The ener gy of this state was 
computed up to the fourth order in the strong coupling expansion in and its mass is 
given by subtracting the ground state energy (jHn|) from it 

nip = Ee- E^ = ela{^ + %a^ -h7%a^). (31) 

The lattice chiral condensate in the staggered fermion formalism may be obtained by consid- 
ering the mass operator M(x) = —l/{Na)Y,x=i{~^Y'^x'^x and evaluating its expectation 
value on the perturbed state > generated by applying to \ip >, >= {i/j > +\p^^^ > 
~^\p^ip^ >. A direct computation of the lattice chiral condensate is given in and yields 



<P^WV^__1/1 32„._^153g„., (32) 

< Pi>\P^ > a\2 



We here compute the ratio between the lattice chiral condensate ()32|) and the lattice pseu- 
doscalar mass ()31|) and equate it to its continuum value 



2 \ \-16z + 768^2 



.ei^V i + 8^- 576^2 27r' 



(33) 



where z = a'^ = A 4, . Eq. ilH^ is true only when Pade approximants are used, since the 
l.h.s. holds only for z ^ 1, while the r.h.s. provides the value of the ratio between the chiral 
condensate and the pseudoscalar mass obtained when z = 00. Using the fourth power of 
eq.(jnSl) in order to construct the [0,2] Pade approximant for the l.h.s. of eq.(|HH|). one gets 



XlV 256 z"" 



(34) 



1 + 384^ + 58368^2 • 

One may now take the limit z —>■ 00 obtaining = l/(228t^) and, equating it to its 

continuum counterpart, one gets an equation for t 

(35) 



228t4 V27r. 

/^From eq. (j35p this one gets t = 1.03 for the lattice light velocity; the result lies 3% above 
the exact value thus showing that the strong coupling evaluation of the chiral condensate 
yields a very good result also for this model. 

In conclusion, our large Nc Hamiltonian approach with staggered fermions evidences 
that the possible ground state of strongly coupled lattice QCD are those of a spin Nc/2 
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antiferromagnetic Ising model; choosing one of the two ground states amounts then to the 
spontaneous breaking of the discrete chiral symmetry corresponding to translations by a 
lattice site. As a consequence a non vanishing chiral condensate is found which is the order 
parameter for the spontaneous breakdown of the chiral symmetry. The bare lattice chiral 
condensate < ipip > is computed in the strong coupling region for a general number of colors 
Nc using the method described in ,1]. The result obtained is intensive, i.e. independent on 
the number of lattice links. Then the 't Hooft limit is taken; the ratio between the chiral 
condensate, obtained in this way, and the third power of the p meson mass computed in the 

n 

't Hooft limit in [1] is used to provide an estimate of the lattice light velocity, t. t turns out 
to be very close to the expected value of 1; this evidences that the large Nc limit is a very 
reliable and relatively simple pathway for the evaluation of the chiral condensate also in the 
strong coupling region. 
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